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Asteroseismic modelling of Procyon A: Preliminary results 

G. Dogan^ * A. Bonanno^ T. R. Bedding^ T. L. Campante^'^ J. Christensen-Dalsgaard^ 
and H. Kjeldsen^ 

^ Department of Physics and Astronomy, Aarhus University, Ny Munkegade, DK-8000, Aarhus C, Denmark 

^ Catania Astrophysical Observatory, Via S.Sofia 78, 95123, Catania, Italy 

^ School of Physics A29, University of Sydney, NSW 2006, Sydney, Australia 

* Centro de Astrofisica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

Received XX, accepted XX 

Key words stars: individual (Procyon A) - stars: oscillations - stars: evolution 

We present our preliminary results of the modelling of the F5 star Procyon A. The frequencies predicted by our models are 
compared with the frequencies extracted through a global fit to the power spectrum obtained by the latest ground-based 
observations, which provides two different mode identification scenarios. 
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1 Introduction 

Procyon A is a member of a binary system with a white 
dwarf companion, Procyon B. It is one of the very bright 
stars to the naked eye, and hence it has been very attrac- 
tive for the observers. The observational constraints which 
we adopted are summarized in Section 2. It has also been 
of asteroseismic interest for a long while (see Arentoft 
et al. 2008, for a summary), with a solar-like power ex- 
cess in its spectrum first reported by Brown et al. (1991). 
However, there has been no agreement on the individ- 
ual oscillation frequencies. Several authors have investi- 
gated the structure and evolution of Procyon A through 
an asteroseismic approach (e.g., Guenther & Demarque 
1993; Barban et al. 1999; Chaboyer et al. 1999; Di Mauro 
& Christensen-Dalsgaard 2001; Eggenberger et al. 2005; 
Provost et al. 2006, Bonanno et al. 2007), but there has 
been a need for more accurate frequencies. Recently, the star 
was observed through a multi-site campaign by eleven tele- 
scopes for more than three weeks (Arentoft et al. 2008). The 
frequency analysis is described by Bedding et al. (2010). 
They presented results from different approaches of fre- 
quency extraction: Iterative sine-wave fitting and global fit- 
ting to the power spectrum. In the former method, a sine 
wave is fitted to each mode one after the other while the 
corresponding sinusoid is subtracted from the time series 
at each step. This is repeated until the signal to noise ra- 
tio of the remaining power is lower than a given threshold. 
This method was used for frequency extraction of ground- 
based radial velocity data before (see, e.g., the analysis on 
solar-like star /3 Hyi by Bedding et al. (2007)). In the latter 
method, the goal is to find an overall fit to the power spec- 
trum for all the frequencies, mode heights, and linewidths 
simultaneously, using some prior knowledge of oscillation 



properties as constraints. A similar implementation of this 
method was previously applied to space-based data (see, 
e.g., frequency analysis of CoRoT star HD 49933 by Beno- 
mar et al. (2009)). In this work, we adopted this Bayesian 
approach which provided us with two mode identification 
scenarios, referred to as Scenario A and B, with different 
posterior odds (for a detailed discussion, see Bedding et 
al. 2010). We chose the output of this analysis in order to 
test both scenarios. Note that Scenario B was favoured by 
Bedding et al. (2010; see that paper for a discussion), but 
here we test both scenarios. 

2 Methods and Tools 

We have adopted the following properties for the position of 
the star in the H-R Diagram: Teg = 6530 ± 90 K (Fuhrmann 
etal. 1997) and log(L/L0)= 0.85 ±0.06 (Steffen 1985). We 
note that there are several revised values for luminosity in 
the Uterature, such as log (i/L© ) = 0. 840 ±0.018 derived by 
Jerzykiewicz & Molenda-Zakowicz (2000) using the Hip- 
parcos parallax and total absolute flux; however, we chose 
to scan a wider range, which largely covers most of the re- 
vised values. A similar argument applies also to the choice 
of the effective temperature. We have not put an additional 
constraint on the radius for the time being, although we do 
compare the stellar mean density infeiTed from the analysis 
with the value 0.172 ± 0.005 Pq, obtained from the mea- 
sured radius using the angular diameter 5.404 ± 0.031 mas 
(Aufdenberg et al. 2005) and the revised Hipparcos parallax 
(284.56 ± 1.26mafQ, van Leeuwen 2007), together with the 
adopted mass 1.463 ± 0.033 Mq, which is the mean of the 
two different astrometric determinations (Girard et al. 2000, 
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' Note that there is an error in Bedding et al. 2010 (Section 9); the value 
they give for the revised pai'allax is actually the original one. We also note 
that the uncertainty on the revised parallax is larger than the original but is 
presumably more reliable. 
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Gatewood and Han 2006). For the metallicity of the star, we 
allowed a wide range covering the 0.05 dcx iron deficiency 
suggested by Allende Prieto et al. (2002), and we used the 
solar mixture from Grevesse & Noels (1993). 

We have computed stellar models with two differ- 
ent evolutionary codes: ASTEC (Aarhus STellar Evolu- 
tion Code) (Christensen-Dalsgaard 2008a) and GARSTEC 
(Garching Stellar Evolution Code) (Weiss & Schlattl 2008). 
The method we used is to compute several grids of stan- 
dard models scanning through a parameter space formed by 
varying the mass, the initial metallicity at the stellar sur- 
face, Zi/X\, and the mixing-length parameter, a, where the 
mixing length is defined as £ = aHp, Hp being the pressure 
scale height. So far, we have searched within the follow- 
ing ranges: M = 1.42 - 1.52 Mq, Z-JX-, = 0.0204 - 0.0245, 
Yi = 0.26 - 0.31, and a = 1.6 - 1.9. Here Xu Yi, and Zu are 
the initial mass fractions of hydrogen, helium, and the ele- 
ments heavier than helium, respectively. We have computed 
the models without taking into account diffusion, convective 
overshooting, or rotation. 

The Aarhus adiabatic pulsation package (ADIPLS) 
(Christensen-Dalsgaard 2008b) has been used to calculate 
the frequencies of the models having properties that are in 
agreement with the observations. We have then compared 
the model frequencies with the observed frequencies. We 
selected the models that minimize the following x^'- 

where f°^^(n), and u™°'^^^{n) are the model, and the ob- 
served, frequencies with spherical degree I and radial order 
n, and a{i'i^^{n)) represents the uncertainties in the ob- 
served frequencies. 

3 Results 

The results from the two different stellar evolution codes 
are similar; hence we present some of the selected models 
computed with ASTEC. The so-called echelle diagrams of 
the best models for both of the scenarios, chosen without 
applying any near-surface corrections are shown in Figs 1 
and 2, with their parameters summarized in Table 1. One 
plots the echelle diagrams using the frequency modulo the 
large frequency separation, Av, in the horizontal axis. In 
order to allow easy comparison with the diagrams shown by 
Bedding et al. (2010), we use the same value, Ai^ = 56 /xHz, 
in our diagrams. 

The behaviour of the frequency differences between the 
models and the observations (shown in Figs 3 and 4) are 
quite different from that in the Sun. Kjeldsen et al. (2008) 
showed that the difference between observed and model fre- 
quencies of the Sun can be fitted by a power law, which can 
also be employed to correct the model frequencies for near- 
surface effects in other solar-like stars, such as (3 Hyi and 
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Fig. 1 Echelle diagram of the selected model for Scenario 
A. Open symbols represent the model frequencies, while the 
filled symbols with the uncertainties correspond to the fre- 
quencies extracted from the observations. Circles, triangles, 
and diamonds are used for the modes with spherical degree 
/= 0, 1, and 2, respectively. The vertical dot-dashed line cor- 
responds to the value of Av. 
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Fig. 2 Echelle diagram of the selected model for Scenario 
B. Symbols are used in the same way as in Fig. 1. 

a Cen A. However, in Procyon, we cannot justify the ap- 
plication of such a surface correction to yield a significant 
improvement in the fit, since the frequency differences do 
not follow the power-law behaviour (see also Figs 16 and 
17 of Bedding etal. 2010). 

4 Discussion and Conclusion 

We can argue, if Scenario A is the correct one, that the pre- 
dictions of the stellar evolutionary models match the ob- 
servations quite well; however, a surface correction for the 
model frequencies seems not to be needed, unlike in the 
solar case. Therefore, the effects of different near-surface 
properties on the high frequencies might be cancelling out 
in Procyon; this deserves further investigation. 

If, on the other hand. Scenario B is correct, there seems 
to be no good agreement between the models and the ob- 
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Table 1 Parameters of the best models 



Parameter 


Scenario A 


Scenario B 


M/Mq 


1.50 


1.50 


Zi/Xi 


0.0235 


0.0245 


Yi 


0.266 


0.290 


Age (Gyr) 


1.83 


1.51 


R/Rq 


2.058 


2.067 


P/Pe 


0.1721 


0.1698 


L/Lq 


6.565 


7.286 


refl(K) 


6446 


6603 


a 


1.8 


1.6 




0.1585 


0.1995 




3.79 


24.46 



* Mass fraction of hydrogen remaining in the centre of the star 
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Fig. 3 The difference between the radial (/ = 0) frequen- 
cies from the observations and the selected model for Sce- 
nario A. The indicated uncertainties are those from the data 
analysis. 

servations, even in the low-frequency region, which means 
that there is something incompatible in the structure of the 
models. 

We have presented preliminary results from our on- 
going work. To come to a conclusion we need to extend 
our analysis. Effects of inclusion of overshooting and use 
of different treatments of convection should be analysed. In 
addition, the suspected mixed mode reported by Bedding et 
al. (2010) could help us distinguish between the two scenar- 
ios, and put further constraints on the age, and the chemical 
composition. 

Although from the point of view of modelling. Scenario 
A seems to be less problematic, we cannot yet strongly 
favour either of the scenarios; therefore, it is difficult to set 
accurate constraints on the stellar properties. Nevertheless, 
either of the cases suggests that Procyon is quite different 
from the Sun, which provides a very good opportunity to 
test our understanding of stellar structure and evolution. 
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Fig. 4 Same as Fig. 3, but for Scenario B 
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